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Abstract 

We give a mathematically rigorous derivation of Ehrenfest's equations for 
the evolution of position and momentum expectation values, under general 
and natural assumptions which include atomic and molecular Hamiltonians 
with Coulomb interactions. 

J* 1 INTRODUCTION 

> . 

A basic result in the physics literature asserts that the mean position and momen- 
tum of a quantum system in M. d with Hamiltonian 

oo 



evolve "classically" 



i> ■ j t ( x i)4>(t) = ^( p j)m> ( 2 ) 
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|tt>«., = (3) 



(Ehrenfest's equations, [Ehr27j ). Here ip(t) 6 L 2 (R d ) denotes the state of the 
quantum system at time t, assumed to evolve via the time-dependent Schrodinger 
^ | equation 

ij t m = Hm, (4) 

(A)^^ stands for the mean (or 'expected value') (ip(t), Aip(t)) of a linear operator 
(or 'observable') A on L 2 (M. d ), (/, g) — f Rd f g is the usual L 2 inner product, X and 
P are the position and momentum operators (see (jHJ), §2$) below), and the potential 
V : R d -> R and the force -W : R d -> R d act by multiplication. The evolution 
of mean values as in @, ([3]), besides being of interest in its own right, plays an 
important role in the study of quantum-classical coupling in molecular dynamics 
(see e.g. |BNS96(lGKZC04j ). Also, a priori estimates related to © are utilized 
in the proofs of important results in scattering theory [SS87, Gra90. |Per93|. IDG97] . 

The heuristic justification, which can be found in any quantum mechanics text- 
book, goes as follows. When H and A are hermitean, formally differentiating the 
mean of A and substituting into the Schrodinger equation yields 

j t {A) m =i{[H,A]) m , (5) 
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where [H, A] denotes the commutator HA—AH. When H is of the form fl} and A is 
a component of the position or momentum operator, formal evaluation of the com- 
mutator gives ([2]), ©• For interesting discussions in the physics literature on the 
shortcomings of the formal argument see [ADVGD011 lADVOOal lADVOObl IHil73j . 
Also, it is straightforward to make Ehrenfest type equations rigorous under suf- 
ficiently stringent assumptions, such as: H is self-adjoint and A is bounded and 
leaves the domain of H invariant; or: H is a Schrodinger operator with smooth po- 
tential, A is a component of the position or momentum operator, and ip(t) belongs 
to Schwartz space for all t, the Schwartz norms being uniformly bounded when t 
belongs to a bounded interval. 

To our knowledge, a rigorous version of j2]), ([3|) under satisfactory assumptions, 
which include in particular the basic atomic and molecular Hamiltonians with 
Coulomb interactions, is so far missing from the literature. It is however readily 
obtained by combining known results (in particular, Hun66> RS78 ) with standard 
functional analytic arguments, as we point out in this note. 

Making rigorous sense of the terms appearing in the abstract Ehrenfest equation 
([5]) when H and A are unbounded self-adjoint operators requires, in particular: 

(1) Well-definedness of the expected value (ip(t), Aip(t)) for all t. When the latter 
is interpreted as an inner product, this means invariance of D(A) under the 
propagator e" ltH of eq. ((4]). 

(2) Differentiability of (ip(t), Aip(t)) with respect to t. 

(3) Well-definedness of the expected value of the commutator HA — AH for all t, 
in some suitable sense. 

As turns out, (1) is essentially sufficient for (2), (3), and the validity of eq. ([5]): 

Theorem 1.1. (Abstract Ehrenfest theorem) Let H and A be two densely defined 
linear operators on a Hilbert space TL such that: 

(HI) H : D{H) —> TL is selfadjoint, A : D(A) — > TL is hermitean 
(H2) e-' ltH leaves D(A) f~l D(H) invariant for allteR 

(H3) For any ip £ D(A) n D(H), sup teJ ||Ae~ ifcff Vo|| < oo for J C K bounded. 
Then for ipo £ D(A) (~l D(H), the expected value (A)M t \, ip(t) := e~ ztH ipo, is con- 
tinuously differ entiable with respect to t and satisfies eq. the right hand side 
being defined as the following quadratic form: 

([H, A])^ := (Hili, Aip) - (Alp, Hip) (ip £ D(A) n D(H)). (6) 

Remarks 1) Assumption (H2) is needed to make the expected value (A)^^ 
well defined (see (1)), and (H3) will be useful in establishing its differentiability 
(see (2) and the proof below). Note that (H2) and (H3) hold automatically when 
A is relatively bounded with respect to H . 

2) Definition ((6]) allows us to make sense of the right hand side of ([5|) without the 
need to consider the (often complicated) domains of the composite operators AH 
and HA, let alone proving their invariance under e~ itH . Note that these domains 
need not be dense when A and H are general self-adjoint operators, and do not 
even contain C§°(M. d ) in basic examples from physics (see below). Of course, when 
ip belongs D(AH) n D(HA) (where D(AH) := {ip £ D(H)\Hip £ D(A)}), © 
reduces to the classical definition (ip, (HA — AH)ip). 

Example Let H be the hydrogen atom Hamiltonian and A the momentum 
operator, H = -\A + V(x), V(x) = -l/\x\, x £ M 3 . OnC£°(IR 3 ) we have i [H, A] = 
-W(x) = -x/\x\ 3 . But -x/\x\ 3 does not map C£°(R 3 ) into L 2 (R 3 ), because the 
singularity at zero behaves like \x\~ 2 and hence |W </>||| 2 = / Rd |VV '(x)(p(x)\ 2 'dx = 



2 



oo for all cj) G C£°(R 3 ) with 0(0) ^ 0. 

3) It is a deeper fact is that if A is in addition assumed to be self-adjoint then (HI) 
and (H2) imply (H3). A proof of this fact, as well as a counterexample showing 
that hermiteanity does not suffice for this implication, will appear elsewhere [FSj . 



4) The main point in the proof is a weak convergence argument which exploits the 
bound in (H3); see Section [2j 

Next, we describe our ensuing results on equations (J3J). Let H be a Hamiltonian 
of the form 

^ 1 B 2 

H = T+V(x), T=-^_^-2, m u ..,m d >0, H = L 2 (R d ), D(H) = H 2 (R d ), 

(7) 

(for conditions on V such that H is well defined and self-adjoint on H 2 (R d ) see 
below) and let A a component of the position or momentum operator, 

{Xjip^x) = Xjip{x), D(Xj) = {ip G L 2 (R d ) | / \ Xj \ 2 \ifj(x)\ 2 dx < oo}, (8) 

Pj = ~, D{Pi) = {if) G L 2 (R d ) \pMp) G L 2 (R d )}, (9) 

(where if) denotes the Fourier transform of if), defined on smooth rapidly decaying 
functions by if>(p) = J Rd e~ w ' x ip(x) dx). Clearly, Xj and Pj are hermitean. In order 
to guarantee well-defmedness and self-adjointness of H on D(H), by the Kato- 
Rellich theorem it suffices to assume that V : R d — > R is a real-valued locally 
integrable function such that 

if) i > Vip is relatively bounded w.r.t.T with relative bound a < 1. (10) 

By this one means that for all if) G iJ 2 (K d ) the function Vip belongs to L 2 (R d ) and 
\\Vip\ \ < a\\Tif)\ \ + C a \\if)\\, for some constants a < 1 and C a . Prototypical are the 
electronic Hamiltonian of a general molecule, d = 3N, x = (yi, ..,yjv), Uj G R , 

and the full electron-nuclei Hamiltonian of such a molecule, d = 3(N + M), x = 
(yi,..,yjv } -Ri,..,-RM), Vi el 3 , R a G R 3 , 

^ = ^-y^A fi + V , ^ , . (11) 

a=l Q l<a</3<M 1 Q pl 

Here the M Q > and Z a > are the masses and charges of the nuclei, and atomic 
units have been used so that h = 1 and the electrons have mass 1 and charge —1. 

Corollary 1.2. (Time evolution of expected position and momentum) 

Let H be given by Q), with V : R d — > R feemg a real-valued locally integrable 

function satisfying the Kato-Rellich condition mU\) . Then: 

(i) (Xj)tPM is continuously differ entiable with respect to t for any ipo G D(Xj) n 
D(H), and satisfies the equation 

j t {Xj) m = i((Hm,x 3 m) - {x 3 m,mt)) 
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(ii) {Pj)^p(t) is continuously differentiable with respect to t for any ipo 6 D(H), 
and satisfies the equation 

j t ( p i)m = i((Hm,PMt)) - (p 3 m,Hm))- 

Proof It suffices to check that the hypotheses of Theorem If .11 are satisfied. 
Clearly, (HI) holds for both Xj and Pj. In the case of Xj, (H2) and (H3) are 
satisfied by the results of |Hun66[ IRS78] . As regards Pj , (H2) and (H3) are satisfied 
since D(H) n D(Pj) = D(H) 1 Pj is relatively bounded with respect to H (i.e. 
\\Pjip\\ < C{\\H^\ \ + \\tp\\) for all ip € D(H)), as is easily deduced from the Kato- 
Rcllich condition, and HV'WII an d II^V'WII are time- invariant. 

Note, however, that so far we have not fully derived the classical Ehrenfest equations 
for position and momentum, since it remains to be verified that the right hand sides 
in (i) and (ii) agree with the classical right hand sides in ((2]), J3}. In particular, in 
order to recover the classical expression {—j^r-)^{t)> i- e - the quadratic form of the 
multiplication operator corresponding to the j th component of the force, additional 
assumptions on the potential are needed. 

Theorem 1.3. Let H be given by where V : R d — > K is real-valued, belongs 
to the Sobolev space W^(R d ) of locally integrable functions with locally integrable 
weak derivatives, and V and yj\ VV| satisfy the Kato-Rellich condition U(J\) . Then: 

(i) (Xj)^,^ is continuously differentiable with respect to t for any ipo £ D(Xj) n 
D(H), and satisfies (0). 

(ii) {Pj)<f){t) is continuously differentiable with respect to t for any ipo G D(H), 
and satisfies (0). 

An important technical point regarding the assumptions of Theorem 11.31 is that 
the Kato-Rellich condition is not required for VV, but only its square root. This is 
related to the fact that VV is only needed as a quadratic form, see Sec. [3l Hence 
the Hamiltonians (fTT|) , (fTTjl satisfy the assumptions of the theorem. This is because 
the square root of the gradient of a typical term in V looks the same as the term 
itself, |V y Jyi— i? Q ,|^ 1 | 1/ ' 2 = \yi — -R Q | _1 . Note that for these Hamiltonians, VV itself 
fails the Kato-Rellich condition (it does not even map Cq° to L 2 , see the example 
in Remark 2) above. 



2 Proof of abstract Ehrenfest theorem 

Proof. We begin by recalling the properties of the strongly continuous one-parameter 
unitary group e~ ttH generated by a self-adjoint operator H (see e.g. (RSSOJ. For 
all t, e~ ltH leaves D(H) invariant and commutes on D(H) with H; moreover for 
any V'o £ D(H), t t— > i/j(t) = e~ ltH ipo is a continuously differentiable map from M 
to H and satisfies i-^ip(t) = Hip(t) for all t. 
In order to show the existence of 

— (tp(t),Aip(t)) =hm , 

dt fc->o h 

we use the decomposition 

W+h), A^(t+h)) - {j,{t), Avp(t)) = / ip{t + h)-jj(t) \ | / ij{t + h)-vp(t) 

h \ h I \ h 

and show the existence of the limits of the two terms on the RHS separately. Since 
^(^(t + h) — ip(t))^ — iH^p{t) strongly in H for h — > 0, the second term converges 
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to i(Hip(t),Aip(t)). Moreover, the first term goes to —i(Aip(t), Hip(t)) provided 
we can show that Aip(t + h)^Aip(t) weakly in Ti for h — > 0. To this end, fix 
t € 1 and choose a sequence {hj} C K. satisfying hj — > 0. Then by (H3), the set 
{Aif}(t+hj)} is bounded in Ti. By weak compactness of the unit ball, there exists a 
subsequence (again labelled by hj) such that Aip(t+hj) — f for some / G Ti. We 
claim that / = Aip(t). To see this, choose any <p G D(A) and calculate using the 
weak convergence of Aip(t + hj), the hermiteanity of A, the continuity of ip(t) in t 
and the fact that ip(t) G D(A) 

(/, <P) = lim (AiP(t+hj), <P) = Um <^(t+/i i ),A0> = (^(t),A0> = (AiP(t),cP). 

hj—fO hj—>0 

Since D(A) is dense in 7i, the claim follows, and since the argument is valid for all 
sequences hj — > 0, it follows that Aip(t + h) — Aip(t). This completes the proof of 
differentiability of (A)^,^ and of eq. ([5]). It remains to show that the derivative 
2lm (Hijj(t), Aip(t)) is continuous in t. Indeed, for h — > we have Hip(t + h) — > 
Hip(t) strongly in Ti and, as just shown, Aip(t+h) — Aip(t) weakly in 7i, completing 
the proof. □ 



3 Derivation of the classical Ehrenfest equations 

Proof of Theorem 11.31 (i) Thanks to Corollary 1 1.21 we only need to evaluate the 
abstract expected value ([6]) of the commutator and show that 

i((Hip, Xjip) - {Xjip, Hip)) = — {ip, Pjip) for all ip G H 2 (R d ) n D{Xj). (12) 
V / rrij 

For functions ip G Co°(K") this follows from an elementary calculation. The general 
case will follow from an approximation argument, but due to the presence of the 
unbounded operator Xj , a little care is needed. 

First, consider functions ip G H 2 (EL d ) with compact support. Clearly it suffices 
to approximate ip by a sequence ip E of C§° functions in such a way that the four 
terms appearing inside the inner products, Hip £ , Xjip s , Pjip £ and ip s , converge in 
L 2 to the corresponding terms for ip. Choose a ball Br(0) containing the support 
of ip. Consider the following standard approximation obtained by mollification: 
^Pe(x) = (Xe * = J Rd Xe(x - y)ip{y)dy, where Xe{x) = e~ d x{x/e) , e G (0, 1), 

X G C °°(R d ), X = outside Bi(0), f Md x = 1- Then (see e.g. [Eva98] l iP £ G C °°, 
■0 e — > ^ in if 2 . Since if, Pj and I are continuous operators from H 2 to L 2 (in case 
of if this follows from (TTUl) ). ^jV'e — * Pjip an( A — * i> m Finally, 

supp^/> £ , supp?/> C S_r+i(0), and is bounded on the subspace of L 2 functions 
with support in Br + i(0), so Xjip e — > Xjip in L 2 . This establishes ([6]) for compactly 
supported if 2 functions. 

Finally let ip be a general function in D(H)C\D(Xj) = H 2 nD(X 3 ). Let x G Cfi° 
with x(0) = 1. Then ipn{x) := x(x/R)ip(x) is a compactly supported if 2 function, 
so (TT2l holds for V'fl by the previous step, and it is straightforward to check that 
ipR —* ip in H 2 . Moreover Xjipn converges to Xjtp in L 2 since XjipR — {Xjip)n 
and X,-i/j G L 2 . Consequently HipR, Xjipn, PjipR and ipR converge in 1? to the 
corresponding terms for ip, establishing (TT2|) in the general case. 

Proof of Theorem ll.3l (ii) As in the proof of (i), it only remains for us to evaluate 
the abstract expected value ([6]) of the commutator and to show that 

i((Hip,Pjip) - (PjiP,Hip)) = (ipi-^ip), (13) 

* v ' S » ' 

=:Q(V>) ~Q*W 
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for all tp G H 2 (R d ). We start by considering ip G C£°(R d ). In this case, an 
elementary calculation shows that 

QW) = [ V(x)^(x)\ 2 dx. (14) 

jR d 3 

Due to the assumption V G W^(R d ) and the fact that \ip(x)\ 2 G C^°(R d ), we can 
integrate by parts to obtain QTSfy for all ip G C£°(R rf ). 

To establish (filjj) for all r/> G -ff 2 , by the density of Cfi° in i? 2 it suffices to show 
that the quadratic forms Q and Q* are continuous on H 2 . As regards Q, this follows 
from the fact that H and Pj are continuous operators from H 2 to L 2 . As regards 
Q*, this follows by writing 

and noting that "0 i— > i— > are continuous operators from H 2 to L 2 by the 

relative boundedness of / with respect to T and the fact that T : H 2 — > L 2 is 
continuous. The proof of Theorem 11.31 is complete. 
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